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SUMMARY 


Piloted simulation studies of cockpit traffic display concepts require the 
development of representative traffic scenarios. With the exception of 
specific aircraft interaction issues, most research questions can be addressed 
using traffic scenarios consisting of prerecorded aircraft movements merged 
together to form a desired traffic pattern. Prerecorded traffic scenarios 
have distinct research advantages, allowing control of traffic encounters with 
repeatability of scenarios between different test subjects. This paper 
describes a technique for generation of prerecorded jet transport traffic 
scenarios suitable for use in piloted simulation studies. Individual flight 
profiles for the aircraft in the scenario are created interactively with a 
computer program designed specifically for this purpose. The profiles are 
then time-correlated and merged into a complete scenario. This technique has 
been used to create traffic scenarios for the Denver, Colorado area with 
operations centered at Stapleton International Airport. Traffic scenarios for 
other areas may also be created using this technique, with appropriate 
modifications made to the navigation fix locations contained in the flight 
profile generation program. 


INTRODUCTION 

For the past several years, NASA Langley Research Center has been involved in 
a joint project with the Federal Aviation Administration (FAA) evaluating the 
Cockpit Display of Traffic Information (CDTI) concept. Langley's role in this 
project has been to study cockpit systems integration issues and pilot roles 
associated with using CDTI in the operation of jet transport aircraft. 
Operations in terminal areas, where traffic congestion is most severe and the 
potential benefits and hazards associated with traffic displays are most 
critical, have been the focus of the work at Langley. The primary research 
technique utilized for this effort has been piloted simulation. 

A critical element in a piloted simulation study involving CDTI is the nature 
of the air traffic scenario utilized in the study. The complexity and realism 
required of a traffic scenario is a direct function of the research objectives 
of a particular study. Part-task experiments, such as those conducted at 
Langley involving in-trail self-spacing with CDTI (ref. 1 and 2), may be 
accomplished with a very minimal air traffic control environment. More 
elaborate experiments, involving full-systems simulation with the test crew 
subjected to a full workload cockpit environment, require a correspondingly 
more complete traffic scenario. This paper will describe a traffic scenario 
generation technique developed for such a full-systems experiment. 

TRAFFIC ENVIRONMENT 

The traffic environment modeled was an area within approximately a 100 n.mi. 
radius of the Denver Stapleton International Airport (see figure 1). 
Technically, all this airspace is the responsibility of the Denver Air Route 
Traffic Control Center ( ARTCC ) • A specific area around the primary airport is 
re-delegated to the Denver Terminal Radar Approach Control (TRACON). The 
ARTCC' s and TRACON' s airspace is further divided, horizontally and vertically, 
into control sectors. Figure 1 illustrates the respective boundaries of the 
high and low altitude sectors of the Denver ARTCC. Included in the figure are 
the locations of several VORTAC (Very high frequency Omni-Range with an 


associated TACtical air navigation transmitter) stations which reside in or 
near the simulated area. Figure 2 is a graphical representation of the 
horizontal boundaries of the Denver TRACON airspace when runway 26L is used 
for instrument approaches. The four metering fixes for Denver are shown in 
figure 2 (KIOWA, KEANN, DRAKO, and BYSON) . The vertical delineation of the 
airspace is shown in figure 3, with the threshold of runway 26L as the zero 
reference in the figure. A table giving a summary of the major airspace 
segments and the significant events which occur in each segment is included in 
figure 3. 

A typical sequence of events for a jet transport aircraft inbound to land at 
Denver Stapleton Airport would begin with air traffic control (ATC) issuing 
instructions for initial descent from cruise altitude when the aircraft is 
approximately 80 to 100 n.mi. from the airport. As the aircraft nears the 
lower vertical boundary of the particular ARTCC sector, coordination is 
initiated for a "handoff," that is, a transfer of control responsibility to 
the next controller. When coordination is completed, the pilot is instructed 
to tune his radio to the frequency of the accepting controller. 

Once the aircraft enters the TRACON airspace (crossing one of the four 
metering fixes), it receives vector headings and altitude assignments which 
merge it into the sequence for the landing runway in use. A final handoff is 
made from the TRACON to the tower controller before the aircraft crosses the 
last navigation fix (outer marker) on the published approach procedure. 

Departure procedures are typically the reverse. Before an aircraft begins to 
taxi, the pilot requests and receives a "clearance" which is an approval of, 
and/or changes to, the requested route of flight and altitude. This is 
followed by taxi instructions to the departure runway in use, takeoff 
clearance, and a transfer of control to the departure radar controller. 

The departure radar controller, usually located in the TRACON, first 
establishes positive radar identification of the aircraft. He then issues 
vector headings and altitude assignments until the aircraft is clear of other 
known traffic and the pilot can assume his own navigation on course. As the 
aircraft nears the boundary of the TRACON' s airspace, it is handed off to the 
adjoining facility. 

Certain conditions were prestated for this particular study. Weather 
conditions were assumed which would require total reliance on the cockpit 
instruments for navigation and aircraft control. Based initially on 
conversations with Denver TRACON personnel, runway 26L was designated the 
landing runway and runways 35L and 35R for departures. Also, a maximum 
traffic density of 35 landings per hour was set. It was further assumed that 
all traffic in the scenarios would be jet transport aircraft. 

After the test scenarios had been designated, traffic data were obtained from 
the Denver TRACON for the months of January through March of 1983. These 
representative data showed that during the sample period instrument weather 
conditions were present 9.9% of the time. During those periods of instrument 
weather, the runway 26L/35 arrival/departure combination was in use 52.6% of 
the time. Also, the mean arrival rate on runway 26L during the Instrument 
Flight Rules ( IFR) periods was 35.9 aircraft per hour, with 89% of the 
aircraft being jet transport types. 
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CREATION OF TRAFFIC SCENARIO 


The creation of a complete traffic scenario is a multistep, iterative 
process. First, the general characteristics of the scenario must be 
defined. Next, individual profiles of aircraft flying the desired approach 
and departure paths in the scenario are created using the computer program 
described in Appendix A. Finally, the individual profiles are time-correlated 
and merged together into the complete scenario. A detailed description of 
this process is given in the following sections. 

Definition Phase 

The first step in creating a traffic scenario is to define the general 
characteristics of the scenario. The major consideration in this definition 
is that the scenario can be designed to provide a specific traffic pattern for 
the test subjects to encounter while flying the simulator through the 
prerecorded scenario. The nature of this traffic pattern is dictated by the 
research objectives of the study which will be using the scenario. The 
desired pattern must conform to the traffic environment being simulated, which 
in this case is present-day Denver under instrument weather conditions. With 
these considerations in mind, the characteristics of the scenario can be 
defined. 

The arrival or departure routing which the simulator will be flying must be 
chosen. The next step is to define the geometry of any specific encounters 
desired between the simulator and any particular traffic aircraft. This may 
be for an in-trail following case, a crossing situation or any desired 
encounter which could occur in the simulated traffic environment. The general 
routing of the encounter aircraft must then be defined from some appropriate 
initial conditions, through the encounter condition on to some final condition 
such as landing or leaving the terminal area. Finally, the flight paths of 
the remaining background traffic must be defined such that the density and 
sequencing of traffic is consistent with the simulated environment. The 
number of aircraft chosen for background traffic must be sufficient to fully 
encompass the simulator within the traffic pattern without having extraneous 
traffic which would burden the real-time simulation. This is not critical at 
this stage since unnecessary traffic will be eliminated during the merge 
phase. 


Profile Generation Phase 

Once the general routings of the simulator and the traffic aircraft have been 
determined, the individual flight profiles of each aircraft must be 
generated. A computer program written specifically for this purpose is 
described in Appendix A. This program provides for computation of the flight 
profile of a jet transport aircraft using interactive or predefined commands, 
which emulate present-day air traffic control techniques. The program was 
designed to enable rapid generation of aircraft flight profiles which are 
realistic for the simulated environment and, at the same time, can be 
controlled to exhibit the desired encounter geometry for the subject research 
study. In order to provide this flexibility, the program was restricted to 
modeling the flight characteristics of the jet transport aircraft without 
specifying any air traffic control strategy. All the commands necessary to 
generate a complete flight profile consistent with the chosen traffic 
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environment were provided. However, the operator of the program must 
determine the proper sequence of commands in order to create a realistic 
profile. Once a profile has been created, a time history of the pertinent 
state variables which define the profile can be saved on a data file for later 
use in the traffic scenario. In addition, a data file of commands used in 
creating the profile is generated. This command data file enables the user to 
recreate the same profile without the need to interactively input all the 
commands again. Fine tuning of the profile can also be accomplished by 
editing the command data file and then rerunning the flight profile program 
with the edited command file. 

The profile generation phase requires an in-depth understanding of the air 
traffic control techniques employed within the simulated traffic 
environment. The first step in creating the profiles is, therefore, to review 
the techniques used by air traffic control for the desired situation. 

Ideally, the assistance of air traffic controllers from the subject facilities 
should be obtained. In addition, the Air Traffic Controller’s Handbook and 
the Airman’s Information Manual (ref. 3 and 4) provide rules and guidelines 
which should be observed. 

Interactive generation of a profile using the computer program described in 
Appendix A is quite straightf orward. The user will be prompted for initial 
conditions for the aircraft. Command options will then be provided allowing 
the user to specify new airspeed or Mach, altitude, heading, landing gear 
position, and ILS approach to runway 26L. Once the desired command options 
(one or more) are selected, the user specifies the segment distance for the 
aircraft to fly before new command options may be entered. At any point 
during the generation of a profile, when the command option menu is displayed, 
the user may review the profile by selecting a plotting option which 
graphically depicts the flight profile up to the current conditions. The user 
may then resume generation of the profile, save the profile data, or start 
over. Throughout this process, the commands input by the user are saved on a 
data file to enable rapid reconstruction of the profile without requiring 
interactive inputs. The actual profile data are only saved at user request 
following the completion of the profile generation. These data consist of 
time, altitude, groundspeed, airspeed, ground positions, heading, and bank 
angle. The data are in a format to be read by the scenario merge phase 
computer program discussed in the next section. 

Merge Phase 

Once the flight profiles of the individual aircraft have been created, they 
need to be merged together into a single complete scenario. This is done by 
time-correlating the flight profiles of the traffic aircraft with the profile 
representing the desired flight path of the simulator in order to produce the 
traffic geometry as specified in the scenario definition phase. Ground 
position time histories of the traffic aircraft are compared with that 
representing the simulator in order to determine a time bias for when the 
traffic profile should begin in order to arrive at the desired location at the 
desired time. The traffic profiles are adjusted by the time biases determined 
in this manner, and are then combined to form a time history of all traffic 
based on the scenario time of the baseline profile representing the 
simulator. In practice, it has been found convenient to use the final 
approach fix (outer marker) of the landing runway as the geographical position 
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for time-correlating the approach traffic. This provides a common point which 
all approach aircraft cross. The time when the simulator crosses the outer 
marker is defined as the baseline scenario time. The desired time for each 
approach aircraft can be referenced to this baseline in order to obtain the 
scenario time of arrival for the aircraft , which provides the desired approach 
spacing and traffic density (landings per hour). The time biases for the 
approach aircraft are then determined by subtracting the actual time they 
would cross the marker from the desired scenario time of arrival for that 
aircraft. Departure and/or conflict traffic must be time-correlated based on 
the desired geometry at some arbitrary geographical location in relation to 
the reference profile of the simulator. A computer program was developed to 
generate the merged scenario, given the individual profiles and time biases. 

A companion program provides plotting capability to view the traffic geometry 
of the merged scenario at any desired time in the scenario. 

Sample Scenario 

Figures 4a through 4d illustrate a traffic scenario created using this traffic 
generation technique. The figures represent "snapshots" of the horizontal 
traffic situation at various times in the scenario. Each figure is centered 
at the threshold of runway 26L in Denver with magnetic north being up and east 
to the right. The map features are the same as those described in the 
Interactive Plotting section of Appendix A. The aircraft representing the 
simulator which is to fly through the scenario (referred to as "reference 
aircraft") is depicted as a circular symbol as identified in figure 4a. The 
traffic aircraft are depicted as squares. The direction of travel for each 
aircraft in the scenario is given by straight line path predictor vectors 
extending from each aircraft symbol. The numbers adjacent to the predictor 
vectors of the traffic aircraft refer to the aircraft identification numbers 
in Table I. The initial condition of the scenario is shown in figure 4a, with 
the reference aircraft crossing the ELBEC intersection at cruise altitude 
enroute to Denver. Figure 4b shows the traffic pattern with the reference 
aircraft just entering the TRACON airspace, approximately 30 n.mi. from the 
Denver Airport. A potential conflict situation between the reference aircraft 
and a departing aircraft is shown in figure 4c. Finally, figure 4d shows the 
traffic situation when the reference aircraft crosses the outer marker on 
final approach to the runway. 

Scenario Definition .- The sample scenario was developed for use in a CDTI 
simulation study in which the test subjects in the simulator would be given 
the task of performing in-trail self-spacing on another aircraft from cruise 
altitude all the way to landing. The traffic environment would be Denver 
under instrument weather conditions. In addition to the self-spacing task, 
another requirement for the scenario was that a different aircraft should pass 
as close to the simulator as possible without causing an actual loss of 
required IFR separation (3 n.mi. horizontal or 1000 ft vertical). This 
encounter should appear to the simulator crew on their CDTI display as a 
possible conflict. Additional background traffic would be needed to include 
all possible aircraft that would either be within range of the CDTI display in 
the simulator or within the same ATC sector receiving radio communications on 
the same frequency as the simulator. Given these requirements, the aircraft 
in the scenario could be defined. 
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The flight profile which the simulator would follow must allow for an extended 
period of in-trail following of another aircraft, as well as the opportunity 
for a close encounter with a third aircraft. Typical operations at Denver 
involve long lines of aircraft crossing one of four entry cornerposts into the 
terminal area. This situation is ideal for in-trail following from any 
direction into Denver. The choice for the actual routing of the simulator 
was, therefore, made to coincide with the best opportunity for an encounter 
with a third aircraft. At Denver, this situation exists when aircraft 
departing the terminal area must pass underneath the airspace of aircraft on 
the downwind leg of the approach to the landing runway. The most common 
runway configuration at Denver is landings on runway 26 and departures on 
runway 35. Given this situation, the most likely conflict will occur between 
approach traffic entering the terminal area from the southwest corner (BYSON 
on figure 2) turning downwind for runway 26L, and departure traffic from 
runway 35L turning 180 degrees left for a south departure (figure 4c). 

Initial conditions for the simulator were, therefore, chosen on a jet route 
southwest of Denver which would bring the aircraft through the BYSON 
cornerpost enroute to Denver. The actual starting point would be the ELBEC 
intersection (see figure 4a) where the simulator would be in cruise just prior 
to initiating its descent into Denver. These conditions would also be used 
for the lead aircraft which the simulator crew would be following with their 
CDTI • The encounter aircraft would be a southbound departure taking off from 
runway 35L. 

The remainder of the traffic in the scenario were defined by sketching what a 
typical landing pattern would look like with the simulator and lead aircraft 
in several positions along their approach. Aircraft would be brought in from 
all approach corridors and merged into a single stream landing on runway 
26L. Starting positions for aircraft approaching Denver from the three 
cornerposts other than BYSON were the holding fixes just outside of the 
cornerpost waypoints. An additional profile from the southwest was also 
defined which merged into the simulator’s route prior to arriving at the 
holding fix for the BYSON approach. In addition to the south departures, an 
east departure was also defined which would be within range of the reference 
aircraft CDTI but would not be a conflict. 

Profile generation .- Having defined the scenario, the individual profiles of 
the aircraft were generated using the program described in Appendix A. Seven 
individual profiles were generated. One profile, starting at ELBEC and ending 
at touchdown on runway 26L, would be used for both the simulator and lead 
aircraft which the simulator would be following. One profile each was created 
for the approach patterns arriving from the other three cornerposts. These, 
and one from the BACCA intersection merging at SHREW and crossing BYSON, were 
generated to represent the remainder of the arrival traffic. Two profiles 
were generated for the south and east bound departures. Air navigation 
charts, the Air Traffic Controllers Handbook, the Airman’s Information Manual, 
and discussions with airline pilots and air traffic controllers from the 
Denver facilities were all utilized to insure that realistic airspeed, 
altitude, and heading commands were used in generating these profiles. 

Merged scenario .- The individual profiles of the aircraft were time-correlated 
and merged into the scenario illustrated in figures 4a through 4d. Table I 
lists the individual aircraft in the scenario with their starting and ending 
positions, the location used for merging the profile, the desired and actual 
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times to this location, and the time bias needed to shift the profile to the 
desired merge time. For example, aircraft number four is the approach traffic 
which the simulator will be following to landing. The profile starts at ELBEC 
intersection and ends on runway 26L. The merge location chosen is the outer 
marker for runway 26L. The desired arrival time at this location is 1230 sec, 
which is 100 sec prior to when the simulator would arrive if it followed the 
exact profile listed as the reference aircraft. Since the actual time that 
aircraft four would arrive at the merge location is 1330 sec, a time bias of 
-100 sec must be applied to correct this arrival time to the desired time. In 
other words, the profile of aircraft four must be advanced 100 sec in order to 
sequence the aircraft into the traffic scenario. The desired arrival times 
for the remainder of the approach traffic were chosen to provide adequate 
spacing between consecutive aircraft and also result in the proper number of 
arrivals per hour as specified in the scenario definition. 

REAL-TIME SIMULATION CONSIDERATIONS 

The final, merged traffic scenario is utilized in the real-time simulation 
program as a sequential input data file which is read at regular time 
intervals during a simulation run. The input data file is created from the 
merged traffic scenario file, with conversions applied to match the simulation 
coordinate reference and time interval at which the real time program will 
read the input file. For example, the basic merged scenario data file is 
referenced to the runway 26L landing threshold with data at two second 
intervals. The DC-9 simulation program required navigation data referenced to 
the Denver VORTAC available at one second intervals. The simulation input 
data file was created by interpolating on the scenario data file to provide 
the necessary one-second update rate, with corrections applied to the data to 
transform the coordinate reference and match that used in the simulation 
program. The resulting input data file was then ready to be used by the 
simulation program without requiring conversions or interpolations while the 
real-time program was executing. 

Two special provisions were incorporated in the DC-9 real-time simulation 
program to enhance the flexibility of the input traffic data file. First, a 
time bias specification for the traffic file was available while the 
simulation program was in a reset mode. This provision allowed advancing the 
traffic file by any desired time interval in order to adjust the initial 
positions of the traffic. Starting location for the simulator could then be 
chosen to be at any position along the designated approach path, with the 
traffic then adjusted to be at the proper position. This was especially 
useful for recreating and demonstrating particular traffic situations. In 
addition, the same traffic file could be used for both approach and departure 
scenarios simply by changing the initial time bias on the traffic. The second 
provision added to the simulation program specified a large table of flight 
identifiers to be assigned to the individual aircraft in the scenario. 
Individual simulation runs were assigned a particular offset value in this 
table used to correlate a flight identifier with an aircraft in the data 
file. This allowed the identical scenario to be used in different data runs 
with different flight identifiers on the traffic without changes to the 
traffic file or tables in the simulation program. This feature was most 
useful in masking an identical traffic scenario by merely changing the call 
signs of the aircraft involved. 
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The ATC coramuni cations with the aircraft in the traffic scenario were 
simulated by creating scripts for the controller to utilize during the "live" 
simulation runs. For each test condition, a traffic data file, time bias, and 
flight identifier offset were established. A printout of the traffic data 
file was used to determine the ATC communications which each aircraft would 
need to receive during the scenario. The time when the communication takes 
place, radio frequency it would occur on, and nature of the communication were 
then recorded. All communications to all aircraft were combined into one 
time-sequenced script for the controller to use during the simulation run. 

Any conflicts in communications were resolved by minor modifications to the 
individual flight profiles involved. During the actual simulation runs, the 
controller was provided with the running scenario time and radio frequency 
which the simulator cockpit had tuned. The controller would then issue or 
respond to any communication which would occur on the frequency tuned by the 
cockpit simulator. By using variations in his voice, the controller was able 
to mimic the responses of the traffic aircraft and the different ATC sector 
controllers during this series of CDTI experiments. An electronic voice 
modifier or several people reading scripts for the traffic could be used to 
provide more realism in the ATC partyline communications. 
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APPENDIX A 

FLIGHT PROFILE COMPUTER PROGRAM DESCRIPTION 
INTRODUCTION 

This appendix describes the computer program used to create the individual 
flight profiles of the aircraft used in generating the traffic scenarios. The 
program was written in FORTRAN 5 for the CDC CYBER 175 computer system at NASA 
Langley Research Center. The program capabilities/ including flight path 
command inputs/ interactive plotting and data file generation are discussed in 
this appendix. A technical description of the equations of motion, 
mathematical models, and control logic used in the program is also provided. 

PROGRAM CAPABILITIES 

This program was written to provide a rapid means of creating realistic flight 
profiles of jet transport aircraft. It provides a set of command inputs 
allowing the user to specify the characteristics of the profile in a step-by- 
step process. The user may enter the command inputs interactively for each 
segment of the profile or from a taped data file. Interactive command inputs 
can be saved on a data file for recreating the profile. An interactive 
plotting capability is provided for viewing the horizontal and vertical 
trajectory parameters of the profile. When the user is satisfied with the 
profile, the trajectory parameters can be saved on a time history data file 
for use in the merged traffic scenario. 

Trajectory Commands 

A flight profile is created as a series of variable-distance segments, each 
defined by a list of trajectory commands. The segments are created in 
sequence with the user specifying the trajectory commands and distance to be 
flown during each successive flight segment. The available commands are 
airspeed, altitude, heading, landing gear, and ILS for 26L. The control logic 
of the program will M fly” the airplane model in such a manner to achieve the 
desired trajectory commands in a realistic fashion. 

Airspeed commands are input as either indicated airspeed or Mach number. 
Altitude commands are entered as feet above mean sea level. Heading commands 
are entered as navigation headings or radar vectors. Navigation headings will 
be flown as a ground track with airplane heading corrected for wind. Radar 
vectors will be flown as airplane heading with ground track allowed to vary 
with the wind. Landing gear selection will toggle the gear up or down 
depending on current status. Selection of ILS for 26L will activate logic in 
the program which will capture and track the localizer and glideslope for 
landing on runway 26L in Denver. This option should not be selected until the 
airplane has been positioned for ILS intercept. 

Trajectory commands may be input to the program interactively or through a 
data file. A new profile is typically generated using interactive inputs. 

The program will prompt the user for all necessary inputs during the 
interactive process. The steps involved in creating a profile in this manner 
are given below. 
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1. Initial conditions of the airplane are specified by the user. 

2. Status of the airplane is displayed. User is prompted for new 
trajectory commands. 

3. User inputs new trajectory commands and specifies a segment distance 
to be flown. Only those trajectory commands which change for the new 
segment need to be entered. 

4. If a segment distance of zero is entered, the user will be given the 
options to plot the trajectory data up to this point on the screen, 
save the data, continue with the profile, start a new profile or exit 
the program. Continuing the profile returns the program to step 2, 
starting a new profile returns to step 1. 

5. Program flies the airplane model for the specified segment 
distance. Control inputs are modified to achieve the trajectory 
commands specified in step 3. The program then returns to step 2. 

The process outlined above provides the user with great flexibility in 
generating flight profiles. At the conclusion of any segment, the profile 
data can be reviewed graphically to insure that everything is proceeding 
properly. The profile generation may then continue, end, or start over at the 
discretion of the user. This allows for correction of errors in a profile 
prior to saving the data, and also permits creation of multiple profiles 
during one session. 

The initial conditions and trajectory commands input to the program during the 
interactive creation of a profile are automatically written to a user- 
specified data file. This data file may then be used in a future session to 
recreate the profile without going through the interactive process. Also 
included in this data file are time and positions of the airplane at each 
segment of the profile. This feature was added to assist in creating radio 
communication scripts with air traffic control instructions for the airplanes 
in the scenarios. The times on the data file where the airplane was given new 
airspeed, altitude, or heading commands were correlated to the scenario times 
in order to determine when such instructions needed to be simulated. 

Another feature of the command data files which proved useful was the ability 
to make minor adjustments in an existing profile. For example, extending the 
downwind leg of a given profile could be done by simply editing the command 
input file and increasing the appropriate segment distances. The profile 
could then be regenerated using the modified command input file. This greatly 
simplified the process of "fine tuning" a given traffic scenario to provide a 
specific traffic pattern. 


Interactive Plotting 

One of the more useful features of this program is the ability to obtain 
interactive plots of the horizontal and vertical trajectory parameters of the 
flight profile. The plotting was accomplished using the PLOT-10 library 
routines contained in the NASA Langley Computer System. A Tektronics graphics 
terminal, or another terminal equipped with PLOT- 10 emulation software, is 
required to utilize this plotting capability. Hard copies of the plots could 
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be obtained from the screen or postprocessed using the plot vector files 
created by the PLOT- 10 routines. 

At the conclusion of any segment during the creation of a profile, the user 
can generate plots of the trajectory parameters. Figure A.1 illustrates an 
example of the plot of the horizontal trajectory of a flight profile. The 
plot shows the ground track of the airplane starting with the initial 
conditions and ending at the final position of the profile. Time marks are 
superimposed on the profile at user-specified intervals as indicated on the 
plot. These time marks are used to correlate the ground position of the 
airplane with time histories of the other trajectory parameters. The 
horizontal map is centered at the threshold of runway 26L in Denver, magnetic 
north being up, with range rings at 10 n.mi. intervals and symbols indicating 
important waypoints and intersections in the Denver area. Also illustrated on 
the map is a dashed line showing the extended centerline for runway 26L. The 
range scale for the plot is specified by the user in terms of n.mi. radius 
from Denver to be displayed. This arrangement is similar to the radar display 
used by air traffic controllers in Denver, and has been found to be most 
useful for providing proper heading vectors to the simulated aircraft during 
the generation of the flight profiles. 

Time history plots are also available for altitude, indicated airspeed, and 
groundspeed of the profile. Figure A. 2 illustrates the altitude and airspeed 
time histories of the profile shown in figure A.1. Altitude or airspeed at 
any point in the profile can be determined by cross-referencing the time marks 
on the horizontal plot with the corresponding time on the plots in figure A. 2. 

Data File Generation 

The end result of this program is a data file containing the trajectory 
parameters of a specific flight profile. The parameters saved include time, 
altitude, groundspeed, indicated airspeed, x position (north) from runway 26L 
threshold, y position (east) from runway 26L threshold, and magnetic 
heading. The data are saved at 2 sec intervals throughout the entire 
profile. Since a complete approach profile can be up to 30 minutes in length, 
this can represent a large amount of data. 

TECHNICAL DESCRIPTION 

The purpose of this program is to generate a time history of the state 
variables describing the flight profile of a jet transport aircraft. This is 
accomplished by integrating simplified trajectory equations of motion over 
various segment lengths along the desired profile. The segments are user- 
defined of variable length, separated by command inputs which specify desired 
characteristics of the profile. These characteristics consist of the 
airspeed, altitude, and heading which the aircraft is to fly. Control logic 
within the program will modify the control variables in the equations of 
motion in order to achieve the specified flight conditions. 

Equations of Motion 

The trajectory equations describing the motion of a point mass over a flat 
earth were used to represent the airplane* These equations ignore all 
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short-term dynamics of the airplane and reduce the math model to a simple drag 
polar and engine thrust representation. These equations are given below: 

mV = Tcos a - D - Wsin y 
(mVcos y) ¥ = (L + Tsin a) sin <|> 

mVy = (L + Tsin a) cos 4> - mg cos y A.1 

. 

h = Vsin y 

. 

x = Vcos y cos Y 
y = Vcos y sin ¥ 


The state variables describing the flight profile are: 

V, Y, Y ' y 

The control variables are: 

T, <(>, a 

where I* and D are functions of a . 

Mathematical Models 

Aerodynamics The use of trajectory equations of motion reduces the 
mathematical description of airplane aerodynamics to simple lift and drag 
relationships. The equations used for these terms are: 

L = C qs 

L A. 2 

D = C qs 
D 


where : 


C_ = C (f) + C (f) a 

Li Li Li 

a o 

Cp. = c: (f) + c (f) c 2 + c 


D D 


0 


Dgear 


(f) 


q = dynamic pressure 
s = wing area 
f = flap setting 
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lift coefficient at zero angle of attack 


slope of lift curve 

drag coefficient at zero lift 

slope of drag coefficient versus square of lift coefficient 


Lift and drag coefficients are solely a function of angle of attack, flap 
setting, and landing gear position. Angle of attack is a control variable 
used by the control logic. Flap setting was programmed as a function of 
indicated airspeed. Landing gear selection was implemented as a main command 
input, and could be raised or lowered as desired by the operator. 

Engine The only engine parameter needed for the equations of motion is 
thrust. Aircraft weight was assumed constant, thus eliminating the need for 
engine fuel flow. Engine dynamics were neglected in this model. Therefore, 
thrust was allowed to instantaneously vary between maximum climb thrust and 
idle thrust. The engine model consisted of tabulated engine thrust for max 
climb and idle power as a function of altitude and Mach number. The data for 
this model were obtained from reference 5 for a JT8D-7 turbofan engine 
operating under standard day conditions. Linear interpolation was used to 
obtain the maximum and minimum thrust available at a given altitude and speed. 

Control Logic 

The control variables available in the equations of motion are thrust, angle 
of attack, and bank angle. Vertical control of altitude and airspeed is 
accomplished by varying thrust and angle of attack while horizontal control 
utilizes bank angle. The control logic implemented in this program was 
designed to provide for changes in flight conditions in a manner consistent 
with typical operation of a jet transport. No attempt was made to model pilot 
response or autopilot logic which might be used to actually fly the aircraft 
through the specified flight conditions. 

Vertical logic .- In order to capture and maintain specified altitude and 
airspeed conditions, four unique vertical flight phases are considered. These 
phases are level flight, climb, descent, and glideslope tracking. The control 
logic examines the altitude error between actual altitude and desired altitude 
in order to determine which flight phase is in effect. The glideslope 
tracking phase is a unique case which is identified by flags passed from the 
ILS tracking routine. Within each flight phase, thrust is set and a desired 
trim flight path angle and trim angle of attack are determined. The commanded 
angle of attack is then calculated using the following equation: 



a 

c 


a t + k 


tY, 


V 


A3 
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where 


a c - commanded angle of attack 

= trim angle of attack for desired conditions 
Y o = previous flight path angle 
Y d = desired flight path angle 
k = gain on gamma error (-1.5 used) 

The desired flight path and thrust settings for the different flight phases 
are given below: 

Level: 

y- 0 

T = thrust required for trim at desired speed and y 
I f more than 3 knots too slow. T = T 

max 

If more than 3 knots too fast. T = T . 

min 

Climb: m m 

T = T 

max 

Y, = trim y at T and desired speed 
d max * 

If more than 2 knots too slow, y. = 0. 

d 


Descent: 


T = 
Y d 


' T . 
min 

= trim y at T . and desired speed 
1 min 

If more than 2 knots too fast, y 


d 


0 . 


Glideslope: 


Y^ - -3 degrees 


T = thrust required for trim at desired speed and y 
I f more than 3 knots too slow, T = T 

max 

If more than 3 knots too fast, T = T 

min 
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Horizontal logic ♦ - Turns are accomplished by adding a bank angle to compensate 
for heading errors. The following equations are used: 


Y = ¥ - ¥ 

err c 

<J> = 5. I ¥ |, limited to 25 degrees 

Y max 1 err * 


where 


Y = heading error 

err 

= commanded heading 

Y = current heading 

<b = max allowable bank angle 

max 


The actual bank angle is changed at a rate of 2 degrees per second whenever a 
heading error exists. Bank angle is never permitted to exceed the maxiumum 
allowable as calculated in the above equation. The sign logic is designed 
such that the aircraft will always turn toward the commanded heading. As a 
result, it is not possible to command the aircraft to turn more than 180 
degrees in either direction with a single commanded heading. If this 
situation is necessary, the turn must be commanded as a series of less than 
180 degree heading changes in the desired direction. 

PROGRAM LISTING 

The following pages contain the program listing for the flight profile 
generation computer program. Although it is written in CDC FORTRAN 5, it can 
be easily modified for use on other computers which support FORTRAN 77 
compilers. The plotting routine utilizes Tektronics Plot-10 subroutines which 
must be linked with the program in order to utilize the plotting capability. 
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PROGRAM F LTPR05 ( INPUT, OUTPUT , TAPE 1 0 , T APE 1 1 , TAPE 12 , T APE 1 3 , 

* TAPE 1 4 , TAPE15, TAPE 16 , TAPE 1 7 , T APE 1 8 , TAPE 1 9 , T APE20 , TAPE21 , 

* TAPE22 , TAPE23, TAPE24 , TAPE25 , T APE26 , TAPE2 7 , TAPE2 8 , TAPE29) 

C 

c ********** ****** *************************** ******* 


PROGRAM FLTPR05 : DATE 08-10-83 


FTN5 VERSION 
ACCEPTS INPUT 
INTERACT I VE 

DATA FILE : 
TAPE 10 THRU 


INTERACTIVELY OR VIA DATA FILE 

: CREATES FILE OF USER INPUTS 
FOR USE AS DATA FILE IN 
SUBSEQUENT SESSION. 

READS USER-SPECIFIED FILE. 

WRITES NEW COMMAND FILE ON SEPARATE TAPE. 

TAPE 1 9 ARE AVA I LABLE 


OUTPUT FLIGHT PROFILE DATA CAN BE SAVED ON 
TAPE20 THRU TAPE29. 


A MAXIMUM OF 5 PROFILES CAN BE CREATED AND 
SAVED DURING ONE SESSION. 

C * * * * * ***************************************** 

c 

C GENERATES APPROACH AND DEPARTURE FLIGHT 
C PROFILES FOR A DC 9 CLASS AIRCRAFT. 

C 

C INCORPORATES LINEAR WIND MODEL 
C 

C VW=VWO+( HP-5330. )*.0012 

C 

C WHERE: VWO'-RUNWAY WIND 
C PSIW-WIND DIRECTION 

C 

C ARE USER INPUTS 

C 

C ************************************************** 

C 

COMMON/TB/TBCL A ( 8 ) ,TBCLA0(8) ,TBCLS(8) ,TBCD0(8) ,TBCDCL(8) , 
*TBCDS(8) , TBCDG( 8 ) , TBF L AP ( 8 ) , XS ( 1000) . YS( 1000) , 

*T I ME ( 1000) , V SAV( 1000) . HPSAV( 1000) , PS I SAV ( 1000) , 

*VCSAV( 1000) 

DIMENSION ENHP ( 10) . ENXMN ( 5 ) ,TTBL(10,5,2),FLAPl(8). 

*F L APH ( 8 ) 

C 

c *X DC-9 DATA ** 

C 

DATA TBF LAP/O. . 1 . , 5 . . 1 5 . , 25 . , 50 . . 50 . . 50 . / 

DATA FLAPl/190. , 165. . 165. . 145. , 135. , 100. , 100. , 100./ 

DATA F L APH/ 500 . ,200. , 180. , 175. . 145. , 135. , 135. , 135./ 
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DATA TBCLA/.096, .1,. 113, .126, .136, .12. .12, .12/ 

DATA TBCLA0/.04 , . 06 , . 19 , . 38 , . 5 , 1 . 32 , 1 . 32 . 1 . 32/ 

DATA T8CLS/0. ,0. ,0. ,0. .0. ,0. ,0. ,0./ 

DATA TBCD0/.021 . .0 305, .0343, .045, .0638, . 141 . .141. . 141/ 

DATA TBCDCL/ .0484 , .045, .0423, .0416, .0409, .0395, .0395, .0395/ 
DATA TBCDG/.0232, .0232, .0219, .0194, .0168, .0105, .0105, .0105/ 
DATA TBCDS/0. ,0.,0.,0.,0.,0.,0.,0./ 

DATA ENHP/0 . , 10000. ,20000. , 30000. ,0.,0.,0.,0.,0.,0./ 

DATA ENXMN/0.0, .2, .4, .6, .8/ 

DATA ( ( ( TTBL ( I , J , K ) , K= 1 , 2 ) , J= 1 , 5 ) , ( = 1 , 4 ) / 

*2 5200 . , 1700. ,22400. ,600. ,20500. ,-200. , 19800. ,-500. , 

*19600. ,-1400. , 

*19800. , 1000. , 17300. , 500. , 16000. ,-100. , 15700. ,-300. , 

*16000. .-900, , 

*14500. , 800. , 12 900. , 300. , 12200. ,0. , 12 100. ,-400. , 

*12400. ,-600. , 

*9600. ,500. ,9100. ,500. ,8800. ,200. ,8850. ,0. , 

*9300. ,-200./ 

S^IOOO. 

C 

C 

c 

c ******** INITIALIZE VARIABLES ********* 

C 

1 ISPOIL-O 
! TRACE-0 
PH I =0 . 

I LS~0 
I LOC— O 
t GL I DE—O 
T I ME ( 1 )— O . 

GAMTOL - .001 
AGA IN = -1.5 
HP I OLE = 15000. 

ROCALT = -3000. 

C 

C ******** INITIAL CONDITIONS ********** 

C 

PRINT 200 

200 FORMAT 1//1X FL I GHT PROFILE PROGRAM — DC 9 ' . 

*/, IX. ' INTERACTIVE = 0, TAPED = 1') 

READ *, I MODE 

C 

IF ( I MODE .EQ. 1 ) THEN 

C 

C<<<< READ INPUT DATA FROM TAPE 2 >>>>>>>>> 

C 

PRINT 601 

601 FORMAT (IX, 'INPUT TAPE NUMBER (10 THRU 19)') 

READ *, I TAPE 

PRINT 602 

602 FORMAT (IX,' OUTPUT TAPE (CANNOT BE SAME AS I NPU I ) ' ) 

READ * , I TAPEW 

REWIND I TAPEW 
REWIND I TAPE 


17 


o o 


c 


600 


READ (I TAPE, 600) 
FORMAT ( A 10 ) 
READ ( I TAPE , *) 
READ ( I TAPE , * ) 
READ ( I TAPE , * ) 
READ ( I TAPE , * ) 
READ ( I TAPE,*) 
READ (I TAPE,*) 
READ ( ITAPE,*) 

IF ( I TRIM .EG. 

READ (ITAPE,*) 

I THR = 2 


T I TLE 

X , Y 

WT, HP , 

I NAV 

VWO, PSIW 
FLAP 
I GEAR 
I TR IM 
1 ) THEN 

ROC 


VC , 


PSIM 


ELSE 

READ (ITAPE,*) I THR 
END I F 
ELSE 


<<<<< INPUT DATA INTERACTIVELY >>>>>>>> 

C 

PRINT 20 3 

203 FORMAT (IX, ’ TAPE NUMBER FOR COMMAND DATA OUTPUT ' , 
* /,* 10 THRU 19 AVAILABLE’) 

READ * , ITAPE1N 
REWIND ITAPEW 
I TAPE— I TAPEW 
C 

PRINT 201 

201 FORMAT (IX, ' INPUT TITLE (10 CHARS MAX)’) 

READ (*,600) TITLE 

PRINT 202 

202 FORMAT (IX, ’INPUT X (NORTH), Y (EAST)’) 

READ * , X , Y 

PRINT 100 

100 FORMAT ( IX, f WT , AL T , KCAS , HEAD I NG ’ ) 

READ * , WT , HP , VC, PSIM 

C 

PRINT 120 

120 FORMAT ( IX, ' 1=NAV HEADING, O^VECTOR’) 

READ *, I NAV 

PRINT 121 

121 FORMAT ( IX, ' W I ND SP EED , D I RECT I ON ’ ) 

READ *,VW0,PS1W 

17 PRINT 101 

101 FORMAT ( IX, ’FLAP SETTING') 

READ * , FLAP 

PRINT 102 

102 FORMAT ( IX, ‘GEAR, O -UP , 1=DN ! ) 

READ *, I GEAR 

PRINT 103 

103 FORMAT ( IX, ’TRIM CONDITION, O^THRUST , 1=ROC ’ ) 

READ * , I TRIM 

IF ( I TR IM . EO. 1 ) THEN 
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104 F0RMAT( IX, 'CL I MB RATE ( FPM) ' ) 

READ * , ROC 

I THR=2 

C 

ELSE 

C 

PRINT 105 

105 FORMAT( IX, 'THRUST (0=IDLE, 1=CLIMB>') 

READ *, I THR 

C 

END I F 

C 

END I F 

C 

C<<< WRITE COMMANDS ON OUTPUT FILE >>> 

C 

' WRITE (ITAPEW.600) TITLE 

WRITE (ITAPEW.*) X, Y. ’ X-NORTH, Y-EAST ' 

WRITE (ITAPEW,*) WT , HP , VC , PS I M , ' WT , ALT , KCAS , HEAD I NG ' 
WRITE (ITAPEW,*) I NAV , ’ I NAV ’ 

WRITE (ITAPEW,*) VWO , PSIW,' WIND SPEED, DIR' 

WRITE (ITAPEW,*) FLAP,’ FLAP SETTING' 

WRITE (ITAPEW,*) I GEAR , ’ GEAR POS (0 = DN,1 = UP)’ 

WRITE (ITAPEW,*) ITRIM,’ ITRIM ( 0 = THRUST, 1 = ROC)' 

IF ( ITRIM . EQ . 1 ) THEN 

WRITE (ITAPEW,*) ROC,' CLIMB RATE (FPM)’ 

ELSE 

WRITE (ITAPEW,*) I THR , ’ THRUST ( 0 = IDLE, 1 = CLIMB)’ 
END I F 

C 

C<<< INITIALIZE VARIABLES BASED ON INPUTS >>> 

C 

X = X * 6076. 

Y = Y * 6076. 

IF ( I NAV . EQ . 1 .OR. VWO .EQ. 0. ) PS I = PSIM 

PSIGC = PSI 
VCOM = VC 
HP COM = HP 
PS I COM •- PSIM 
C 

DO 15 I = 8 , 1 , - 1 

IF ( FLAP .EQ. TBFLAP ( I ) ) GO TO 16 

15 CONTINUE 

16 I FLAP = I 

CALL FLAPS ( F L AP L . F L APH , VC , I F L AP ) 

C 

THETA--- ( 1— . 0000068 7 *HP ) 

DEL=THETA**5.256 

XMN=5*( ( ( ( 1 .+. 2* (VC/661 . 5 ) **2 ) **3 . 5- 1 . J/DEL+1 . )**. 286-1 . ) 
XMN=SQRT(XMN) 

03=1481 . *DEL*XMN*XMN*S 

VT=2 9 . Q4*SQRT (518. 7*THETA )*XMN 

I F ( VWO . EO.O. )VG=VT 

CALL WIND (HP, VWO, VT, VG.PS I .PSIW. PSIM, I NAV ) 

PRINT 90, VT , VG, PSIM, PSI 
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90 FORMAT ( IX . ' VT-’ . F 10 . 2 . ’ VG= ’ , F 10 . 2 , ' . PS I M= ' , F 10 . 2 . ' PSI=’. 

*F 10.2) 

I F ( INAV.EQ. 1 )PSICOM=PSIM 

CALL ENG ( HP . XMN . ENHP . ENXMN , TTBL , TMAX , TM I N ) 

I F ( I THR. EQ.O)T=TMIN 
I F ( ITHR.EO. 1)T=TMAX 
C 

C ** INITIAL AERO TERMS ** 

C 

CLAO=TBCLAO( I FLAP ) 

CLA=TBCLA( I FLAP ) 

CDO=TBCDO( I FLAP) 

I F ( I GEAR . EO . 1 )CDO=CDO+TBCDG( I FLAP) 

CDCL=T8CDCL ( I FLAP ) 

C 

C ** INITIAL TRIM ** 

C 

IF< ITRIM.EQ. 1 ) GAMMA— AS I N ( (R0C/(VT*101 .28) ) ) 

CALL TR I M( I TR I M , T , WT , OS , CL AO , CL A , CDO . CDCL , GAMMA .ALPHA , 

* PH I ) 

ROC=VT*SIN(GAMMA)*101 .28 
GAMMA=GAMMA*5 7 . 3*VT/VG 

C 

C 

C ** I N I T I AL I ZE MAI N LOOP ** 

C 

I CNT=2 

XS(2)=X/6076. 

YS ( 2 )=Y/6076 . 

VSAV( 2 )=VG 
VCSAV ( 2 ) -VC 
HPSAV(2)=HP 
PS I SAV ( 2 )=PS I M 
T I ME ( 2 ) =0 . 

RANGE = O. 

C 

C*********: MAIN LOOP xr*********** 

C 

2 CFLAP=TBFLAP( I FLAP) 

C 

IF ( IMODE .EQ. 1 ) THEN 

C 

C<<<<<< READ TIME HEADER WITH DUMMY ARGUMENT >>>>>> 

C 

READ (I TAPE. 600) DUMM 
READ (I TAPE, 600) DUMM 
READ (I TAPE, 600) DUMM 
END I F 
C 

C<<<<<< WRITE TIME HEADER >>)>>>> 

C 

WRITE ( ITAPEW,*) ' ' 

WRITE (ITAPEW,*) ’ *** TIME = ’ , T I ME ( I CNT ) . ’ ***' 

* . ' RANGE = ’ .RANGE 

WRITE (ITAPEW,*) ’ NORTH = XS(ICNT). ' EAST = YS(ICNT) 
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c: 

CALL COMMAND ( HP , VC , PS I M , HPCOM , VCOM , PS I COM , SD I ST , IGFAR. 

*CFLAP, I TRACE , ILS, 1 NAV , VT , VWO , PS I W , PS I GC , I MODE , I TAPE, ITAPEW) 

I F ( SD I ST . EO . O . ) GO TO 5 
DCHECK=0. 

4 I SAV E=0 

3 CALL ENG( HP , XMN , ENHP , ENXMN , TTBL , TMAX , TM I N ) 

C 

IF ( I LS . EQ. 1 ) THEN 

C 

CALL I LS26L ( VT , VG , X , Y , HP , HPCOM, PS I , PS I ft, PS i COM, I LOG , 

* I GL I DE , VftO ) 

C 

END I F 

C 

C 

IF (VCOM . LT. 1. ) THEN 

C 

VCCOM = 

$ 148 1 . *SQRT ( ( DEL* ( ( VCOM**2/5 . +1 . ) **3 . 5-1 . )+1 . )**. 286-1 . ) 

XMNTRM=VCOM 

C 

ELSE 

C 

V C COM= V COM 
XMNTRM = 

$ SORT ( 5* ( ( ( ( 1 . + . 2*( VCCOM/661 . 5)**2 )**3. 6-1 . )/DEL+1 . )**. 286-1 . ) ) 

C 

END I F 
C 

OSTR I M= 1481. *DEl*XMNTRM**2*S 

C 

DH - HP - HPCOM 
DV =s VC — VCCOM 
CLAO = TBCL AO ( I FLAP) 

CLA = TBCL A ( I FLAP) 

CDO = T BCDO ( I FLAP) 

IF (I GEAR .EO. 1) CDO = CDO + TBCDG ( I FLAP ) 

CDCL = T BCDCL ( I FLAP) 

C 

c MA I N CONTROL LOGIC 

C 

IF ( ILS .EO. 1 .AND. I GL I DE .EO. 1 ) THEN 

C 

C ON GL I DESLOPE 

C 

I TR IM = 1 

OGAM = -3./S7.3 * VG /VT 

CALL TRIM( 1 , T , ftT , OSTR I M , CLAO , CL A , COO , CDCL , DGAM, DALP , PH I ) 

IF ( DV .LT. -3. ) T=TMAX 

IF (DV .GT. 3. ) T-TMIN 

C 

ELSE IF ( DH .GT. 50. ) THEN 

C 

C DESCENT 
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I TR J M 


0 


T = TM I N 

IF ( HP .GT. HP IDLE .AND. ROC .LT. ROCALT ) THEN 
I TRIM = 1 

DGAM = ASIN ( ROCALT / (VT * 101 . 28 ) ) 


C 

C 

C 


C 

C 

C 

C 

C 

C 

C 


C 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 


c 

c- 

c 

c 

c 

c 


ELSE I F ( DV .QE. 2 . ) THEN 

LEVEL OFF IF TOO FAST 

I TRIM = 1 
DGAM = 0. 

END DESCENT 

END I F 

ELSE IF ( DH .LT. -50. ) THEN 

CL I MB 

T ^ TMAX 
I TRIM = 0 

IF ( DV .LT. -2. ) THEN 

LEVEL OFF IF TOO SLOW 

DGAM - 0. 

I TRIM = 1 

END CL I MB 

END I F 


ELSE 


LEVEL 

I TR IM - 1 

DGAM = -DH * .001 

CALL TR I M( 1 , T , WT . OSTR I M , CL AO . CL A , CDO , CDCL , DGAM, DAL P , PHI ) 
IF (DV .LT. -3. ) T = TMAX 
I F (DV .GT. 3. ) T = TMIN 

END I F 


75 CONT I NUE 

— FIND TRIM ALPHA AND GAMMA — 

PHI = PH I /5 7 . 3 

77 XCL = ( WT*COS ( DGAM ) /COS ( PH I ) )/QSTRlM 


DALP = ( XCL— CL AO ) /CL A 
t F ( I TRIM . EO. 1 )GO TO 76 
XCD = CDO + CDCL*XCL*XCL 
XDD = XCD*QSTRIM 

DGAM1 = AS I N ( ( T*COS ( DAL P/5 7 . 3 ) — XDD ) /WT ) 

I F ( ABS( DGAM— DGAM1 ) .LT. GAMTOL ) GO TO 76 
OGAM = DGAM1 
GO TO 77 

76 DGAM = DGAM* 5 7 . 3 
PHI = PHI *57 . 3 

— MODIFY ALPHA BASED ON GAMMA ERROR — 

ALPHA = DALP + AGA I N* ( GAMMA-DGAM) 

I F ( ALPHA .GT. 1 2 . ) AL PHA=1 2 . 

I F ( ALPHA .LT. 0.)ALPHA=0. 

IF ( I NAV .NE. 1 ) GO TO 85 

CALL WIND (HP, VWO.VT.VGDUM.PS I GC.PSIW.PS I COM, 1) 

86 CALL BANK( PS I M . PS I COM . PH I ) 

CALL AERO( I FLAP , I GEAR , I SPO I L . ALPHA . CL , CD ) 

X L— CL *QS 
D=CD*QS 

CALL SOLVE ( WT , ALPHA . GAMMA , PS I ,PHI , T , XL , D . X , Y , HP , OR , ROC , 
*VG , VT , PS I M) 

CALL ATMOS( HP , VT , VC , XMN , OS , OEL ) 

CALL W I NO ( HP , VWO , VT , VG , PS I ,PSIW,PSIM, I NAV ) 

QS=QS*S 

CALL FLAPS( FLAPL , FLAPH , VC , I FLAP ) 

I F ( I F L AP . EO . 0 ) GO TO 32 
I SAVE— I SAVE+1 
DCHECK=DCHECK+DR 
RANGE = RANGE + DR 
t F ( I SAVE. LT . 10) GO TO 3 
I CNT = I CNT+ 1 

T I ME ( I CNT ) — T I ME ( ICNT-1 )+2 . 

I F ( I TRACE . NE . 1 )GO TO 87 
PRINT 91 , T IME( I CNT ) 

PRINT 91 , VC, HP . ALPHA , T 
PRINT 91 , PS I , PH I , ROC 
91 FORMAT( IX. 'TRACE 1 ,4F8.2) 

87 CONTINUE 

XS( I CNT )=X/6076 . 

YS ( ICNT)=Y/6076. 

VSAV ( I CNT )=VG 
VCSAV ( I CNT ) =VC 
HPSAV ( I CNT ) =HP 
PS I SAV ( I CNT )— PS I M 
I F ( DCHECK . LT . SD I ST ) GO TO 4 
GO TO 2 
5 CONTINUE 

PRINT 301 
READ *, IPLOTF 

IF ( IPLOTF .EO. 1 ) CALL PLOT (ICNT) 



o o 


c 

PRINT 300 
READ *, I SAV 
IF (I SAV .EG. 1) THEN 
C 

<<<<<<<<< SAVE FLIGHT PROFILE DATA >>>>>>>>>>> 

(OUT = I TAPE + 10 
WRITE ( I OUT , 600 ) TITLE 
DO 25 I =2,1 CN r 

WRITE ( I OUT , * ) TIME(I), HPSAV(I), VSAV<I), PSISAV(I), 
* XS( I ) , YS ( I ) , VCSAV ( I ) 

25 CONTINUE 

PR INT 401 , lOUr 
END I F 
C 

PRINT 400 
READ *, ICONT 

IF ( ICONT . EQ , 1 ) GO TO 2 

300 FORMAT ( IX, 'SAVE DATA ( O-NO , I^YES ) ' ) 

301 FORMAT( IX , ’PLOT DATA (1=YES) f ) 

400 FORMAT ( IX, 'CONTINUE (1=YES) f ) 

401 FORMAT ( IX . 'DATA SAVED ON TAPE', 12) 

32 PRINT 31 

31 FORMAT( IX . ' FLAP ERROR') 

30 PRINT 500 

500 FORMAT( IX , 'NEW CASE (1=YES)*) 

READ * , I NEW 

I F ( I NEW. EO. 1 )G0 TO 1 

END 

SUBROUT I NE TR I M( I TR I M , T , WT , OS , CL AO , CL A , CDO , CDCL . 

*GAM, ALPHA, PHI ) 

C 

Q XXXXXXXXXXXXKXXXXX.X.XXXXXX.XXXX.X.XX.XX.XXX'X.XX.XX.XX.XX.X 

C 

C THIS ROUTINE CALCULATES TRIM ANGLE OF ATTACK 

C FOR : 

C 

C I TRIM = 0 ; FIXED THRUST, VARIABLE GAMMA 

C 

C 1TR1M = 1 ; FIXED GAMMA, VARIABLE THRUST 

C 

C 

ACAL. —0 . 
r 

IF ( I TRIM .EO. 1 ) THEN 
TIN IT = 1000. 

C 

ELSE 

T INI T = T 
C 

END I F 
C 

11 AJNIT-ACAL 
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20 XL=WT-TINIT*AINIT/57.3 
CL=XL/QS 

Cl=Cl/COS( PH 1/57.3) 

ACAL=( Cl— CLAO ) /CL A 
ATEST=ABS ( A I N I T-ACAL ) 

I F ( A TEST . GT . 0 . 05 ) GO TO 11 

CD=CDO+CDCL*CL*CL 

D=CD*QS 

IF( ITRIM.EO.O) GO TO 30 
TC=D+WT*GAM 

TEST=ABS( T I N I T-TC ) /T I N I T 
I F (TEST. LT. 0.005) GO TO 21 
T INI T-TC 
GO TO 11 

21 T=TC 

GO TO 40 
30 GAM=( T— D )/WT 
40 ALPHA-ACAI. 

RETURN 

END 

SUBROUTINE COMMAND ( HP . VC , PS I . HPCOM , VCCOM , HCOM , SD I ST , I GEAR , 
*CF LAP , I TRACE , I LS . I NAV , VT . VWO , PS I W . PS I GC , I MODE. I TAPE , I TAPEW) 
C 

Q * X X**XX*.XXX*X.*X**.X.X*t* X** * XXX XX XXXXOK* **)K*)lf XX XXX XX 

C 

C THIS ROUTINE ACCEPTS INPUT TO ALTER THE 
C FOLLOWING COMMAND PARAMETERS: 

C 

C HP COM (DESIRED ALTITUDE) 

C VCCOM (DESIRED AIRSPEED) 

C HCOM (DESIRED HEADING) 

C TRACE (DEBUG TRACE PRINT OUT) 

C ILS (ILS MODE SELECT) 

C I GEAR (LANDING GEAR SELECT) 

C 

c * XX* XXX XXX* ***.* tt*********************** * 

c 

DIMENSION I COM ( 6 ) 

PRINT 10 , HP . VC , PS I 

10 FORMAT (/. IX, *ALT=' .F8.0./, IX, ' KCAS- ' ,F7.0./.1X.'PSI-',F8.1) 
PRINT 200.CFLAP 

200 FORMAT ( IX FLAPS ‘.F3.0) 

11 CONTINUE; 

C 

IF (I MODE .EQ. 0 ) 

* PRINT 20 

20 FORMAT (/, IX , ’NEW COMMAND ’,/, 1 X 1=KCAS . 2=AL T , 3=PS I ’./ . 

*1X . ' 4— TRACE , 5= I LS26L . 6=GEAR ' ) 

C 

C 

IF ( IMODE .EO. 1 ) THEN 

READ ( I TAPE. *) NC 
ELSE 

READ *.NC 
END I F 
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WRITE (ITAPEW.*) NC . * NEW COMMAND 


C 

IF (NC .EQ. 6) THEN 
C 

I F ( I GEAR . EQ. O) THEN 
I GEAR=1 
PRINT 100 

100 FORMAT ( IX . ’GEAR DOWN * ) 

C 

El SE 

I GEAR=0 
PRINT 101 

101 FORMAT ( 1 X . ‘ GEAR UP*) 

C 

END I F 
C 

ELSE IF (NC .EO. 5) THEN 
IF ( IMODE . EO. 0 ) THEN 
PRINT 81 

81 FORMAT ( 1 X I LS OPTION ( 1=V ES , 0=NO ) ’ ) 

READ * , I LS 
ELSE 

READ ( I TAPE , * ) I LS 
END I F 

WRITE ( I T APEW , * ) ILS,* I LS MODE SELECT . 1 = ON’ 

C 

ELSE IF (NC .EO. 4) THEN 
IF ( IMODE .EO. 0 ) THEN 
PRINT 71 

71 FORMAT ( IX TRACE OPTION (1=YES)*) 

READ * , I TRACE 
END I F 
C 

ELSE IF (NC .EQ. 3) THEN 
IF ( IMODE .EO. 1 ) THEN 

READ (I TAPE.*) HCOM 
READ (I TAPE.*) I NAV 
ELSE 

PRINT 21 

21 FORMAT ( IX , ’NEW HEADING’) 

READ * , HCOM 
PRINT 300 

300 FORMAT (IX.’ 1— NAV HEADING. O^VECTOR’) 

READ *, 1 NAV 
END I F 

WRITE (ITAPEW.*) HCOM. ’ NEW HEADING* 

WRITE (ITAPEW,*) I NAV , * 0 = VECTOR. 1 = NAV* 

C 

IF (I NAV . EO. 1 ) THEN 
PS I GC-HCOM 

CALL W I ND ( HP , VWO . V T . VGDUM , PS l GC . PS I W . HCOM . 1 ) 

C 


C 


END I F 

ELSE IF (NC .EO. 1) THEN 
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IF (I MODE . EQ. 1 ) THEN 

READ ( I TAPE , * ) VCCOM 
ELSE 

PRINT 22 

22 FORMAT ( IX, 'NEW AIRSPEED OR MACH') 

READ * , VCCOM 

END I F 

WRITE (ITAPEW,*) VCCOM, * NEW SPEED' 

C 

ELSE IF (NC .EO. 2) THEN 
IF (IMODE .EO. 1 ) THEN 

READ (I TAPE,*) HPCOM 
ELSE 

PRINT 23 

23 FORMAT ( IX / NEW ALTITUDE' > 

READ * , HPCOM 

END I F 

WRITE (ITAPEW,*) HPCOM,' NEW ALTITUDE' 

C 

ELSE IF (NC .EO. 0) THEN 
IF ( IMODE .EO. 1 ) THEN 
READ ( I T APE , * ) SD I ST 
ELSE 

PRINT 61 

61 FORMAT ( / , ’ SEGMENT DISTANCE') 

READ *, SO I SI 
END I F 

WRITE (ITAPEW,*) SD I ST , ’ SEGMENT DISTANCE' 
RETURN 
C 

END I F 

C 

UU IU 11 

c 

END 

SUBROUT I NE ATMOS ( HP , VT , VC , XMN , Q, DEL ) 

THETAS ( 1 .-. 0000068 7*HP ) 

SIG-THETA**4.2S6 

DEL=THETA**5.256 

RHO=. 00237 7*S 10 

Q=. 5*RHO*( ( VT*1 . 688 )**2) 

XMN=SQRT(Q/( 1481 . *DEL ) ) 

QC=Q* ( XMN* XMN/ 4 . + 1 . ) 

VC-1481 . *SQRT ( ( 1 . *f QC/ 2116. )**. 286-1 . ) 

RETURN 

END 

SUBROUT J NE ENG ( HP , XMN , ENHP , ENXMN , T TBL , TMAX , TM I N ) 

C 

C ************************************ 

C 

C THIS ROUTINE CALCULATES r 
C 

C TMAX (MAXIMUM CLIMB THRUST) 

C 

C TMIN ( IDLE THRUST) 
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c 

c 

C FOR GIVEN ALTITUDE AND MACH NUMBER 
C 

C ******************************* ****** 

C 

DIMENSION ENHP ( 8 ) ,ENXMN(8) , TTBL (10,5,2) 

IHP=2 

10 i F ( HP . L E . ENHP ( I HP ) ) GO TO 11 
I HP- I HP+ 1 

l F ( I HP . GT . 3 ) GO TO 11 
GO TO 10 

1 1 I XMN=2 

20 I F ( XMN . L E . ENXMN ( I XMN ) ) GO TO 21 
I XMN= I XMN+ 1 

I F( I XMN. GT. 4) GO TO 21 
GO TO 20 

21 SLXMN=(XMN-ENXMN( I XMN-1 ) )/( ENXMN ( I XMN )— ENXMN ( I XMN-1 ) ) 
SLHP=( HP— ENHP ( I HP— 1 ) )/(ENHP( IHP)-ENHP( IHP-1 ) ) 

TMAX 1=SL XMN* ( TTBL ( I HP-1 , I XMN , 1 ) -TTBL ( I HP-1 , I XMN-1 ,1)) 
TMI N1=SLXMN*( TTBL < I HP-1 , I XMN , 2 ) — TTBL ( I HP-1 , I XMN-1 , 2 ) ) 
TMAX 2=SL XMN* ( TTBL ( I HP , ! XMN , 1)-TTBL( IMP, I XMN-1, 1) ) 

TM I N2=SLXMN*( TTBL ( I HP , I XMN , 2 ) — TTBL ( I HP , I XMN-1 ,2) ) 

TMAX 1=TMAX 1+TTBL ( I HP-1 , I XMN-1 , 1 ) 

TMIN1=TMIN1+TTBL( IHP-1, I XMN-1, 2) 

TMAX2=TMAX2+TTBL( I HP , I XMN-1 , 1 ) 

TM I N2=TM I N2+TTBL ( I HP , I XMN-1 ,2) 

TMAX=TMAX 1+SL HP* ( TMAX2— TMAX 1 ) 

TM I N=TM I N 1+SLHP* ( TM I N2-TM INI) 

RETURN 

END 

SUBROUTINE BANK ( PS I , PS I COM , PH I ) 

C 

C ********************* **************** 

c 

C THIS ROUTINE COMPUTES A BANK ANGLE 
C TO CORRECT HEAD I NG ERROR. 

C 

C PS I = CURRENT HEADING OR COARSE 
C 

C PS I COM = DESIRED HEADING OR COARSE 

PHI INCREMENTED 1 DEG PER ITERATION 
C 

C PH I MAX LIMITS BANK BASED ON PS I ERR 
C 

Q ********************* *************** 

C 

DPH I =1 . 

PSCOMT=PS I COM 
PS I T=PSI 

I F ( PS I . GT . 180 . )PSI T^PSI T-360. 

I F ( PSCOMT . GT . 1 80 . ) PSCOMT=PSCOMT-360 . 

PS 1 ERR=PSCOMT-PS I T 
I F ( ABS ( PS I ERR ) .LT. .1 )GO TO 15 
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PHIMAX = 5 . *ABS ( PS I ERR ) 

IF ( PHIMAX .GT. 25. ) PHIMAX = 25. 

IF(PSIERR.LT.O. ) GO TO 11 
I F ( PS I ERR . GT . 180. ) GO TO 10 
GO TO 20 

11 I F(PSI ERR. LT.-180. )GO TO 20 
GO TO 10 
15 PS I “PS I COM 
PH 1=0. 

GO TO 30 
10 PHI=PHI-DPHI 

I F ( PH I . LT . ~~PH I MAX ) PH I =~PH I MAX 
GO TO 30 
20 PH I =PH I +DPH I 

I F ( PH I . GT . PH I MAX ) PH I =PH I MAX 
30 RETURN 
END 

SUBROUTINE AERO ( I FLAP, I GEAR , I SPO I L , ALPHA , CL , CD ) 

C 

C **** *********** **** ******************** ****** 

C 

C THIS ROUTINE COMPUTES CL AND CD FOR 
C SPECIFIED FLAP , GEAR AND SPOILER 
C CONF IGURAT ION. 

C 

C {FLAP = INDEX CORRESPONDING TO FLAP 
C SETTING IN TBFLAP ARRAY. 

C 

C I GEAR = O (GEAR UP), 1 (GEAR DN) 

C 

C I SPO I L = O (SPOILERS RETRACTED ) , 1 (EXTENDED) 

C 

c ***** XX** XXXXXXXXXXX****************** XXXXxX* 

C 

COMMON/ T 8/ T BC L A ( 8 ) , TBCLAOl 8 ) , TBCLS<8) , TBCD018) ,TBCDCL(8) , 
* T BOOS ( 8 ) , TBCDG ( 8 ) , TBFLAP ( 8 ) ,XS(1000) , YS( 1000) , 

*T I ME ( 1000) , VSA V< 1000) ,HPSAV( 1000) ,PSISAV( 1000) , 

* VCSAV ( 1000) 

CL=TBCL AO ( I FLAP)+ALPHA*TBCLA( I FLAP ) 

CD=TBCDO( I FLAP )-FTBCDCL ( I FLAP)*CL*CL 
IF( I GEAR . EO . 1 )CD=CD+TBCDG( I FLAP ) 

I F ( I SPO I L . EQ . O ) GO TO 20 
CL=CL FTBCL S ( I FLAP) 

CD=CD+TBCDS ( I FLAP) 

20 RETURN 
END 

SUBROUTINE FLAPS( FLAPL , FLAPH , VC, I FLAP) 

C 

C *****xx****xxxxxxxx.xxx**xxxxxx*xxxxxxx 
C 

C THIS ROUTINE DETERMINES PROPER FLAP 
C SETTING BASED ON NOMINAL AIRSPEED 
C SCHEDULE. 


*************** ************ *********** 
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c 

DIMENSION FLAPLC8) ,FLAPH(8) 

10 VL=F LAPL ( I FLAP ) 

I F ( VC . GT . VL )GO TO 11 
I FL AP= I FLAP+1 
I F ( I FLAP . GT . 8 ) GO TO 20 
GO TO 10 

11 VH=FLAPH( I FLAP ) 

IFTVC.LT. VH)GO TO 12 
I FLAP=I FLAP-1 

I F( I FLAP . GE . 1 ) GO TO 1 1 
20 I FLAP=0 

12 RETURN 
END 

SUBROUTINE SOL VE ( WT , ALPHA . GAMMA , PS I .PHI ,T,XL, 

*D . X . Y , HP , DR , ROC . VG , V T . PS I M ) 

C 

Q ** ***** ***************************** 

c 

C THIS ROUTINE SOLVES THE TRAJECTORY 
C EOUATIONS OF MOTION. 

C 

C ******* ****************** *********** 

90 FORMAT ( IX, 'SOLVE ’.4F10.3) 

DT=. 2 

A=AL PHA/57 . 3 
G=GAMMA/57 . 3 
PS=PS 1/57 . 3 
PH=PH 1/57.3 
XM=WT/32 . 2 
VT=VT*1 . 688 
VG=VG* 1 . 688 

CALCULATE "DOT" TERMS 

VDOT=( T*C0S ( A ) -D— WT*S I N ( G ) ) /XM 
PSDOT=(XL+T*SIN(A) )*SlN(PH)/(XM*VG*COS(G) ) 

GOOT=( XL+T*SI N( A ) ) *COS ( PH )/( XM*VG ) -32 . 2*COS(G)/VG 
XDOT=VG*COS ( G ) *COS ( PS ) 

YDOT=VG*COS ( G ) *S I N( PS ) 

HDOT=VG*SIN(G) 

INTEGRATE 

VT— ( V T + VOOT*DT ) / 1 .688 
VG=(VG+VD0T*DT)/1 .688 
PS I M=PS I M+PSDOT*D T*5 7 . 3 
PS I — ( PS+PSDOT *DT ) *5 7 . 3 
GAMMA=( G+GDOT*DT >*57.3 
DX=XDOT*DT 
DY=YDOT*DT 

DR=SQRT(DX*DX+DY*DY )/6076 . 

X=X+OX 
Y-Y+DY 


30 



HP=HP+HDOT*DT 
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c 

C CORRECT HEADING 

C 

I F ( PS I . GT . 360 . ) PS I =PS I —360 . 

IE (PS! . LT . 0 • ) PS l =PS I +360 . 

I F(PSIM.GT. 360. ) PS I M=PS I M-360 . 

I F ( PS I Mt LT • 0 • ) PS I M=PS I M+360 . 

91 FORMAT (IX,* SOLVE2 * , 4F 10 . 3 ) 

ROC=HDOT*60. 

RETURN 

ENO 

SUBROUT I NE PLOT ( I CNT ) 

C 

c ************** *tt******Xtt*X**ttX**X*tt****X>< 

C 

C THIS ROUTINE USES PLOT 10 FOR INTERACTIVE 
C PLOTTING OF . 

C 

C GROUND TRACK (XLOC) 

C 

C TIME HISTORIES 

C 

C ALTITUDE (HP) 

C SPEED (V) 

C GROUNDSPEED ( VG ) 

C KCAS (VC) 

C 

C ************** * **** xx *x* xxxxxxxx********* 

C 

COMMON/TB/TBCLA ( 8 ) ,TBCLA0(8) ,TBCLS(8) . TBCD0(8) . TBCDCL(8) , 
*TBCDS(8) , TBCDG( 8 ) . TBFLAP ( 8 ) . XS( 1000 ) , YS( 1000) . 

*T I ME ( 1000) , VSAV ( 1000) , HPSAV( 1000) . PS I SAV( 1000) , 

*VCSAV ( 1000) 

DIMENSION X ( 200 ) , Y ( 200 ) 

1 PRINT 100 

100 FORMAT ( 1 X , ' X VAR I ABL E ' . / , 1 X . ' 1=X LOC , 2=T I ME ' ) 

READ *, I X 

I F ( I X . EO . 1 ) GO TO 11 
PRINT 101 

101 FORMAT (IX. ' Y VAR I ABLE ‘ IX , ’ 1=V , 2=ALT ' ) 

READ * . I Y 

11 CALL I N I TT ( 120) 

CALL B I N I T r 

CALL SL I MX ( 300, 900) 

CALL XFRM( 3 ) 

CALL XMFRM( 3 ) 

CALL YFRM( 3 ) 

CALL YMFRM( 3 ) 

I F( I X.GT. 1 )GO TO 10 
XS( 1 )= I CNT— 1 
YS ( 1 )= I CNT— 1 
CALL ZLINE(O) 

PRINT 103 

103 FORMAT ( IX, 'RANGE SCALE (NM)') 
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READ * , SCALE 
XLOW=— SCALE 
XH I GH=SCALE 
YLOW=— SCALE 
YH I GH=SCALE 
PRINT 104 

104 FORMAT ( IX , ’ T I ME INTERVAL FOR TIC MARKS') 
READ *, I T 1C 
I T I C= I T I C/2 

CALL DL I MX ( XLOW , XH I GH ) 

CALL DL I MY ( YLOW , YH I GH ) 

CALL CHECK ( YS , XS ) 

CALL DSPLAY ( YS , XS ) 

I T=1 

I TT=I T 1C 
DO 60 1=2. I CNT 
| TT=I TT+1 

I F ( ITT.LT. ITIOGO TO 60 

I T= I T+1 

X( I T )=XS( I ) 

Y ( I T)=YS( I ) 

I TT=0 

60 CONTINUE 
X< 1 )=l T-1 
Y ( 1 ) = I T-1 
CALL ZLINE(-I) 

CALL SYMBL ( 8 ) 

CALL CHECK ( Y , X ) 

CALL CP LOT ( Y , X ) 

C 

C ** VOR SITES ** 

C 

C DEN, IOC, GUC 
C 

X( 1 )=3. 

Y ( 1 )=3. 

X ( 2 ) =— 1 . 79 
Y ( 2 )=2 . 40 
X ( 3 > =2 7.69 
Y ( 3 ) =— 13.59 
X ( 4 ) =—82 . 10 
Y ( 4 ) =— 9 9 . 13 
CALL SYMBL ( 1 ) 

CALL CHECK ( X , Y ) 

CALL CPLOT(X.Y) 

C 

C ** WAYPOINTS ** 

C 

C ENROUTE: 

C 

C ELBEC. ACREE, SHREW, BENAM, SMITY, RAMAH 
C 

X( 1 )=15. 

Y ( 1 ) = 1 5 . 

X ( 2 ) =— 6 1.24 


32 



Y ( 2 ) =—78 . 99 
X<3)=— 42.54 
Y ( 3 ) =—60 . 93 
X ( 4 ) =-30 . 62 
Y ( 4 ) =— 42.34 
X ( 5 ) =—36 . 14 
Y ( 5 ) =34 . 4 7 
X ( 6 )=39 . 97 
Y ( 6 )=42 . 97 
X ( 7 ) =45 . 9 1 
Y(7)=— 26.68 
C 

C TERMINAL! 

C 

C BYSON, DRAKO, KEANN. TR02E , JASIN, PLOTS. WIFES 
C ALTUR, WAT K I 
C 

X(8)=— 21 .43 
Y<8>=— 28.09 
X ( 9 )=— 25 . 74 
Y(9>=24. 74 
X(10)--22.54 
Y ( 1 0 ) =26 . 04 
X( 1 1 ) =— 9 . 71 
Y( 1 1 )=-9. 99 
X( 12) =— 12.74 
Y ( 1 2 )=12 . 64 
X( 1 3 )= 10 . 26 
Y ( 1 3 )=14 . 10 
X( 14 ) — 1 3 . 14 
Y ( 14) =— 5.70 
X( 15)=5. 35 
Y( 15)=1 . 16 
X( 1 6 ) = 1 3 . 4 7 
Y ( 16 )=3 . 16 
CALL SYMBL ( 3 ) 

CALL CHECK (X.Y) 

CALL CPLOT(X.Y) 

CALL ZLINE(I) 

CALL SYMBL ( 0 ) 
r 

C X* RANGE rings ** 

C 

I FOCAL E .GT. 40.) GO TO 53 
X( 1 )=180. 

Y ( 1 )= 1 80 . 

R=0 . 

DO 50 1=1 , 3 
R=R+ 1 0 . 

TH=0 . 

DO 52 J=2.181 
TH=TH+ .0349 
X(J)=R*COS(TH) 

Y(.J)=R*SIN(TH) 

52 CONTINUE 
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CALL CHECK <X,Y) 

CALL CP LOT ( X , Y ) 

50 CONTINUE 
53 CONTINUE 

** EXTENDED CENTER-LINE ** 

X( 1 )=2 . 

Y ( 1 ) —2 . 

X ( 2 ) =0 . 

Y ( 2 ) =0 . 

X ( 3 ) =1 9 . 5 
Y ( 3 ) —4 . 5 
CALL CHECK ( X , Y ) 

CALL CPLOT ( X , Y ) 

CALL MOVABS( 350 , 700 ) 

CALL TSEND 
ITIC=2*ITIC 
PRINT 105, ITIC 

105 FORMAT ( 1 X , ' T I ME I NT ERVAL^ * # I 5 
GO TO 900 
10 T I ME ( 1 )=ICNT -1 

I F ( I Y . GT . 1 ) GO TO 20 
PRINT 250 

250 FORMAT ( IX, f 1— GNDSPEED , O— CAS ' ) 
READ * . I SP 

PRINT 251 

251 FORMAT ( 1 X , ' YM * N , YMAX , f ) 

READ *,YMIN,YMAX 

CALL DL I MY ( YM I N . YMAX ) 

I F ( I SP . EO. 1 ) GO TO 12 

VCSAV ( 1 ) = I CNT— 1 

CALL CHECK( TIME VCSAV) 

CALL DSPLAY(T I ME , VCSAV) 

GO TO 900 
12 VSA V ( 1 ) ” I CNT™ 1 

CALL CHECK ( T I ME , VSAV ) 

CALL DSPLAY(T I ME , VSAV) 

GO TO 900 
20 HPSAV ( 1 ) =■ I CNT — 1 
PRINT 252 

252 FORMAT ( IX, ! YM IN, YMAX, f ) 

READ * , YM I N , YMAX 

CALL DL IMY(YMIN, YMAX) 

CALL CHECK( T I ME , HPSAV ) 

CALL DSPLAY(T I ME , HPSAV) 

900 CONTINUE 

CALL MOVABS(0 , 700) 

CALL TSEND 
PRINT 102 

102 FORMAT ( IX, NEIAI PLOT ( 1-YES)*) 
READ *, 1C 
I F ( 1C. EO . 1 ) GO TO 1 
RETURN 
END 


* SEC * ) 


r 
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SUBROUTINE I L S26L ( VT . VG . X . Y . HP . HPCOM . PS I . PS I W . PS I COM . I LOC , 
* I GL IDE.VWO.PSIM) 

C 

C ********************* ***** *************** 

C 

C THIS ROUTINE HANDLES I LS APPROACH TO 
C RUNWAY 26L AT DENVER. 

C 

C COMMANDED HEADING (PS I COM) MODIFIED 
C TO CAPTURE RUNWAY EXTENDED CENTERLINE. 

C 

C COMMANDED ALTITUDE (HPCOM) MODIFIED 
C TO FLY GLIOESLOPE. 

C 

C **************************************** 

c 

PS I RNY=257 . 

X THRS=0 . 

YTHRS=0. 

DXTHRS-XTHRS— Y/6076 . 

DXTHRS=ABS( DXTHRS ) 

DYLOC=ABS( DXTHRS* . 23 1+YTHRS— X/6076 . ) 

I F ( I LOC. EO. 1 )GO TO 10 
DEL PS I =AB$( PS I - PS I RNY ) 

I F (DELPSI .GT . 180. ) DELPS I =DEL PS I -180 . 

PS I DOT =5 7. 3*32.2*TAN(25./57. 3)/(VT*1 .688) 

T I ME-DELPS I / ( PS I DOT *3600 . ) 
0YTURN=SIN(DELPSt/114.6)*VT*TIME 
C PRINT 1 10 , DYTURN , DYLOC 

110 FORMAT (IX,' D YTURN= ’ , F 10 . 4 . ' DY LOC= ' . F 10 . 4 ) 

IF (DYLOC. GT.DYTURN)GO TO 20 
I L0C=1 

10 CONTINUE 

IF ( VWO .NE. 0. ) GO TO 1 1 

PS I COM = PS I RNY 
GO TO 12 

11 CALL W I ND( HP, VWO, VT.VG.PS I RNY. PS I W, PS I COM, 1) 

12 IF (DYLOC .GT. .017) GO TO 20 
PSI = PSIRNY 

I F ( VWO . EO. 0. ) PSIM = PSI 

C 

C ** GL I DESLOPE CHECK ** 

C 

20 XTD=ABS(DXTHRS)+1000./6076. 

HPGSL P=53 30 . + .0523*XTD*6076. 

DHP— HPGSLP-HP 

I F( ABS(DHP) . LT . 50. )G0 TO 30 
IGL I DE=0 

I F ( DHP . L T . 0 . ) HPCOM— HPGSL P 
GO TO 40 
30 IGLIDE-1 
HPCOM^HP 

40 I F ( HP . LE . 5330 . ) HPCOM-- 53 30 . 

RETURN 

END 
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SUBROUTINE IN I NO < HP . VWO , VT , VG , PS I , PSIiAI,PSlM, INAV) 

C 

q ********************* ************************* 

c 

C THIS ROUTINE CORRECTS GROUNDSPEED ANO HEADING 
FOR WIND EFFECTS. 

PS I - GROUND COURSE 
PSIM = MAGNETIC HEADING 

INAV = 0 SIGNIFIES HEADING MODE .• 

C GROUND COURSE COMPUTED BASED ON REFERENCE 
C MAG HEADING. 

C 

C INAV “ 1 SIGNIFIES NAVIGATION MODE : 

C 

C MAGNETIC HEADING COMPUTED BASED ON REFERENCE 
C GROUND COURSE . 

C 

C ********************************************** 

C 

I F < VWO . EO . 0 . ) GO TO 40 
VW=VWO+( HP-5330 . )*.0012 
PS I WV=PS I W— 180 . 

I F ( PS I WV . LT , 0 . ) PS I WV=PS I WV + 360 . 

I F ( INAV . EO. 1 )G0 TO 20 
10 DPS I W— ( PS I WV— PS I M) /5 7 . 3 

DPS I =AT AN ( VW*S I N ( DPS I W ) / ( V T + VW*COS ( DPS I W ) ) ) 

PS I —PS I M+DPS I *57 . 3 
GO TO 30 

20 DPS I =AS I N ( VW*S I N( (PSIWV-PSI )/57.3)/VT) 

PS I M=^PS I —DPS 1*5 7 . 3 
DPS I W— ( PS I WV- PS I M ) /57 . 3 
30 VG-V T *COS ( DPS I >+VW*COS ( DPS I W -DPS I ) 

RETURN 
40 VG ■ VT 

PS I =PS I M 

RETURN 

END 

— EOR — 

— EOI/TOP — 

? ? 
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All times expressed in secon 
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Range, n. mi . 



A1 rspace 
segment 

Altitude; 

Typical distance from airport; 
Responsible controller(s) 

I 

Above 24 000 ft; 
Greater than 80 n.mi.; 
High altitude enroute 

II 

20 000 to 24 000 ft; 
Greater than 35 n.mi .; 
Low altitude enroute 

III 

Surface to 20 000 ft; 

10 to 35 n.mi.; 

Approach and departure sectors 

IV* 

7 000 to 10 000 ft; 

5 to 20 n.mi . east; 
Final approach control 

V 

Surface to 6 500 ft; 
0 to 5 n.mi . ; 

Local control 


Significant 

events 


Enroute cruise, 
descent clearance 


Profile descent, 
enroute climb 
clearance 


Vectored to final 
approach zone or 
departure course 


Vectored to ILS 


Landing/Takeoff 

clearance 


•Note: Runway 26L arrival configuration 


Figure 3.- Vertical view of simulated airspace 
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